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I. INTRODUCTION 



The problem of motion of a fluid in its own gravitational field (collapse) refers to key problems of astrophysics. Until 
the end of the 1990's, studies of models of compact objects and a collapse process have been basically concentrated on 
<^ , a consideration of models in which a fluid possesses a positive (or zero) pressure. However, the situation has changed 
Qh' after the discovery of the accelerated expansion of the present universe [H, [2l ■ It is now believed that such acceleration 
Q ] is caused by the presence in the universe of some antigravitating matter with the negative pressure p < — 1/3 e, where 
$_i ' e is the energy density (for a review, see e.g. 0, It is assumed that such matter, which is called dark energy, is 
c/3 ' distributed uniformly in space and provides more than 70% of a total mass of the universe. However, a possibility is 
not excluded that dark energy could affect the formation of the structure of the universe on small scales of the order 
of galaxies, and may be even smaller (for a more detailed discussion of this question, see Ref. [sj). 
. A more intriguing possibility is a description of dark energy and dark matter within the framework of some unified 
J> ' theory by using an effective hydrodynamical description of a dark fluid. One of variants of such a theory is a model 
] of the Chaplygin gas @ . It is assumed in this model that the dark fluid is being described by the polytropic equation 
of state pch ~ — Pch' where p and p arc the pressure and the energy density, respectively, and 7 is a negative constant. 
In a cosmological aspect, such equation of state gives the following scenario of the evolution of the universe: in the 
early stages of the expansion of the universe pch ~ a~'^ (here a is the scale factor) that corresponds to the epoch of 
^^O ■ domination of non-relativistic (dustlike) dark matter. On the other hand, at the present time and later on such fluid 
, acts as a cosmological constant, pch ^ const. This attractive feature allows to use such model for a description of the 
evolution of the universe from the epoch of formation of the large scale structure until the present time. 

Possessing the properties of dark matter, the Chaplygin gas can form compact astrophysical objects such as stars 
and galaxies. Such objects had been already considered earlier. In particular, in Ref. 01 the model of the dark energy 
star supported by the Chaplygin gas was considered within the framework of Newtonian gravity. Choosing various 
^^j ■ values of a central density of such configuration, sizes and masses of dark stars in the early universe were estimated. 
d I In general relativity, such stars were considered in detail in Refs. [1, Q. 

In this paper we consider the motion of the Chaplygin gas in its own gravitational field within the framework of 
Newtonian gravity. Such model can be used for a description of evolution of the dark fiuid inside galaxies. Since the 
sizes of such objects are significantly smaller than cosmological scales then one can use the Newtonian approximation. 
In the paper Q the parameters of the Chaplygin dark star were estimated, taking into account the cosmological 
evolution of the Chaplygin gas. Here we will estimate parameters of a motion of the Chaplygin gas creating a 
compact astrophysical object. To do this, we will consider a self-similar motion of the gas with the above equation 
of state. Such problems for self-gravitating gaseous spheres have been already considered by a number of authors 
(see e.g. earlier works [Tol - [T2j . and also more recent papers [H, [13] and references therein). In these works the main 
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attention was called to the models of collapsing stars with the parameter 7 > 0. It will be shown below that in the 
case of 7 < (the Chaplygin gas) self-similar solutions considered here describe only expanding configurations. 

II. EQUATIONS AND SOLUTIONS 
A. Similarity equations 

We will consider the spherically symmetric motion of the Chaplygin gas in Newtonian gravity. The corresponding 
hydrodynamical equations can be written in the form [l2j | 

— +4Trr^pu = 0, (1) 
dM , 

-gf - P> (2) 
du du _ Idp GM 
dt dr p dr ' 

where u is the radial velocity of the Chaplygin gas, M is the total mass of matter contained in the interval (0,r), p 
and p are the density and the pressure of matter, respectively. In the case under consideration, the pressure p is small 
compared with the energy density p. Equations ([I])-© are invariant under the time reversal operation T(^: 

t — > — t, u — > — M, p ^ p, p ^ p, A/ — > M. 

Thus one can consider only the range of < i < 00. Let us choose the polytropic equation of state of the generalized 
Chaplygin gas in the form 

p = -Kp\ (4) 

where K is some positive arbitrary constant and —1 < 7 < 0. The case 7 = — 1 corresponds to the usual Chaplygin 
gas d). 

In this paper, we will seek similarity solutions of the system of equations (HJ-© with the equation of state (|3]). For 
this purpose, we introduce the following dimensionless similarity variables (for details, see [12l|): 

I, jL3/2j-3n-2 

where k is some dimension constant. Using these similarity variables, Eqs. ([l])-(l3]) can be rewritten in the form 

to(x) = ax^[nx~v)^ (6) 

da dv X V 

inx ~ v)— ct-r ~ ^2 a, (7) 

dx dx X 

^nda dv , . nx — V 

_^Q,7 {^nx-v}— = -(n-l)w-- -a. (8) 

dx dx in — 1 

The similarity variables ([5]) in general suppose that 

K = k{ATiGy-^t^^''+^-^\ (9) 

If we restrict ourselves to the case of constant K (as it is usually assumed at a consideration of the dark Chaplygin 
fluid), the parameter n must be restricted by the additional condition 

71 = 2-7, (10) 

or, taking into account 2 < n < 3. For convenience of performing the mathematical analysis, let us rewrite Eqs. 
([7]) and ([8]) in the following form 



da a 



dx [nx — v)^ + 'ya'* ^ 

dv _ 1 

dx {nx — vy + ja'^~^ 



nx V ix v){nx — v) 

(ri-l)fH a-2> ' 



371 — 2 X 



(nx ~ v) X — V 1 

(71 — l)(7ia; — v)v H a + 2 7 a ' 

371 — 2 X 



(11) 
(12) 
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B. Limiting behavior 

In general, the above equations correspond to a non-autonomous system which should be investigated numerically. 
However, let us estimate first the behavior of the system ((TT|) - ((T^ near the origin of coordinates. Wc will follow Rcf. 
p^ : we expand a and ?; in a Taylor series in the neighborhood of x = 0: 



a{x) — + aix + a2x'^ 



v{x) = Wo + vix + V2x'^ + 



(13) 



Substituting these expressions in Eqs. pT|) and (HH, one can obtain the following limiting solutions for a; <C 1 
parameterized by a* : 



v(x) 
a{x) 



2 al~'' 
—X H 

3 157 

2-7 



n \ X 

3/ 



a 



67 V"*~3'^ 



Note that Eqs. ([7]) and ([8]) have a particular solution with the constant density a 

2 2 2 / 2' 



"=3"' "=3' '"=3r"3'" 



(14) 
(15) 

(16) 



corresponding to the homogeneous distribution of matter and describing expansion in Newtonian cosmology. 

Using (fT4 )) ■ (fTS ]) as initial conditions, let us solve Eqs. ([TT|) . (fT2|) numerically. The results are presented in Fig. [TJ 
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FIG. 1: The graphs of velocities v(x) (the solid lines) and densities a{x) (the dashed lines) at different initial values: a* 
0.66, 0.6, 0.5, 0.1, from top to bottom, for both sets of the curves. 



The numerical analysis indicates two regions of the solutions. In the first region a* < 2/3. In this case there exists 
the critical point in which the numerators and denominators of Eqs. (jlip . (jl2p vanish simultaneously. The position 
of this point can be defined from the conditions: 



-70; 



(7-l)/2 



(nx — v) = 

nx ~ V (x — v)(nx — v) 

(n - l)v + a - 2-i '-^ '- = 0. 

^ ^371-2 x 



(17) 
(18) 
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An asymptotic form of solutions of Eqs. pT|) . (|12p at 2; ^ 1 is defined from the form of the numerical solutions. 
One can see from the numerical calculations that v changes more slowly than x, and the function a'^~^ changes more 
slowly than x'^. It allows to find the following asymptotic solutions: from Eq. (fTTj) we have 



da 2a 
dx nx 



^a;-'/". (19) 



Using this solution, Eq. (fT2|) can be rewritten as 



dv ^ {n-l)v ^ _4_^_2/„ ^ 27^^_i^_2(^_i+„)/„ 
dx nx 3n — 2 



with solution 

Zn - 2 [2(n + 7) - 3] n 

Here A, B are integration constants. Taking into account the relation n = 2 — 7, one can see from the last equation 
that a dominating term is the term with B (if, of course, one omits the case -6 = 0). Then we have 

V = (20) 



Taking into account the restrictions on values of — 1 < 7 < 0, 2 < ?i < 3, one can see from Eqs. (|T9l) and ((20|) that 
the density of the Chaplygin gas is a slowly decreasing function at a; 3> 1: a ^ x""^^"^ at 7 = —1 and a ^ x^^ at 
7 — > —0. At the same time, the velocity v is an increasing function: v ^ x^^^ at 7 = —1 and v ^ x^l"^ at 7 — > —0. 

In the second region of the solutions, when a, > 2/3, there also exists the critical point which is defined by equations 
([T7)) and However, in this case, somewhere on the solution plane of the system of equations ([TT|) -([T ^ . there 

also exists a point where the denominator of Eq. ([TT|) goes to zero but the numerator is not equal to zero. I.e. the 
solution becomes singular, and the function a, as opposed to the region a* < 2/3, is not a decreasing function but 
an increasing one that is also physically unacceptable. 

Because of the above, let us use the regular similarity solutions ([T4| . p5|). (fTQ]) and ([20|) in the region a* < 2/3 for 
purposes of estimation of the physical quantities ([Sj rewritten via the variables (r, i) . Using the definition for x from 
([S]), one can see that the limiting case x 3> 1 corresponds to the initial instant t — > (with the arbitrary finite values 
of r). Then, using the expressions (|19p . ((20)) . we have from ([5]) 

u(r, t) = Bk'/^^\^"-'^/\ p{r, t) = ^^r-'^". (21) 

47rG 

I.e. at the initial instant there is the static inhomogcneous distribution of the Chaplygin gas. On the other hand, the 
limiting case a; ^ 1 corresponds to i — > 00 at some finite r. The corresponding velocity and density of matter from 
© will then be 

I.e. asymptotically, as t — >■ cxj, the model describes the homogeneous distribution of matter which does not depend 
on the equation of state (i.e. on n). If one follows the motion of some layer of matter then after transition to the 
Lagrangian variables, equation (|22p for the velocity u can be easily integrated giving the following expression 



ro 



to 



where tq, to SiVC integration constants describing characteristic parameters of the configuration. In cosmology language, 
this equation corresponds to the motion of dust matter. I.e. within the framework of the model under consideration, 
the Chaplygin gas, creating the compact astrophysical object, describes asymptotically, as t — > 00, a cloud of the 
dark fluid with the homogeneous distribution of matter with the dustlike equation of state. In this case the initially 
inhomogcneous distribution of matter (PT|) smooths asymptotically as t — 0. 
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III. CHAPLYGIN GAS IN GALAXIES 



Since the solutions ([2T|) are self-similar, i.e. they are similar to each other at different times, then one can use 
these solutions for a description of the present distribution of the dark fluid in the universe. Since the Chaplygin gas 
may mimic dark matter, we want to apply the solution obtained above as a dynamical model of dark matter inside a 
spiral galaxy. We assume that dark matter inside a galaxy is a dynamical object evolving in time. We will apply the 
asymptotic solution (PT|) for a description of the dark matter distribution in a galaxy ion x ^ 1 (that corresponds to 
t = const., r ^ our model a dark matter halo evolves in time, and we try to describe such evolution at 

some instant. As an example, let us consider a possibility in the use of the solutions (PT|) for modeling dark matter in 
spiral galaxies. For such galaxies there is the so-called Universal Rotation Curve defining the velocities of the motion 



of stars in a galaxy at any radius. The form of such curve can be represented by the following expression [15 1 
VuRcii) = V{Ropt) 



0.72+0.44 log + 1.6e-"-4(^/^-)- 



Lj (£2 + 0.782)1-43 • + (L/L,) 



0.4 



1/2 

km s"\ (23) 



where x = r / Ropt, Ropt is the radius encompas sing 83% of the total integrated light, L is the galaxy luminosity, and 
logL^/i© ~ 10.4. As was pointed out in Ref. [1^, Eq. (|23p predicts rotation velocities at any radius with a typical 
accuracy of 4%. The first term in this formula gives the contribution from a stellar disc, and the second one ~ the 
contribution from a dark halo. For our purposes, let us use the second term. Then we have from ([23|) the following 
distribution of velocities which is defined by the dark halo 



VB'Mr/R^pt) ^ V{Ropt) Jl.6e~0AiL/L,)_^^^^_^ km s'^. (24) 

y + 2.25 (L/L^) 

Let us now estimate velocities of test particles moving in the field of gravity of the dark halo created by the 
Chaplygin gas. The corresponding profile of the velocity distribution will then be given by the formula 



y-r(.) = ^^. 

The expression for the mass distribution of the dark fluid follows from Eq. ([6]). Taking into account Eqs. (O and 
(PT|) . the expression for the mass profile at x 1 will be 

M{r)^ /'^fe'^" ^(3n-2)/n 

^ ' {3n - 2)G 

and the corresponding velocity distribution is 



For convenience of comparison of this expression with formula (|24p , let us rewrite it via the variables x and Ropt as 

r)D("-l)/" _ -, 

vSi:r'ii) = yiRopt)Di"', where D= °g . , m = (25) 

y \n.opt ) n 

Further, choosing the parameters D and m, one can compare Eq. (|25p with formula (|24p . As an example, let us take 
L/L* = 0.75. The corresponding parameters will then be: D w 0.59 and m w 0.67 (see Fig. [2]). Since the parameter 
D contains the integration constant A then its value can be chosen arbitrarily. The value m sa 0.67 corresponds to 
n « 3, i.e. 7 w — 1 (see eq. (fTO)) ') that corresponds to the usual Chaplygin gas. 

Summarizing the obtained results, we have considered the self-similar motion of the Chaplygin gas with the equation 
of state (H]). For values of the parameters used in creating the cosmological models with the Chaplygin gas (— 1 < 
7 < 0), it was shown that only expanding solutions do exist. This is caused by the nature of the Chaplygin gas 
having a negative pressure. Within the framework of the model under consideration, the numerical analysis was 
performed, and also the analytical estimations of the behavior of the solutions were made for the limiting cases when 
the similarity variable a; ^ 1 and a; <C 1 (see Eqs. (|2ip and ([2^ . respectively). The case of a; ^ 1, corresponding 
to the static distribution of the Chaplygin gas, has been applied for the evaluation of the velocity profile of stars in 
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FIG. 2: Comparison of observational data 
parameters D ~ 0.59 and m ~ 0.67. 
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(the dots) with the curve H25|l for the Chaplygin gas (the sohd line) with the 



galaxies. Comparing the observational curve for a dark halo and the curve for the Chaplygin gas it was 

shown that they are in good agreement with each other at the corresponding choice of the fitting parameters D and 
m (see Fig. [5]). 
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